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The interplay between between gauge- field winding numbers, 6- vacua, and the Dirac operator 
spectrum in hnite-volume gauge theories is reconsidered. To assess the weight of each topological 
sector, we compare the mass-dependent chiral condensate in gauge field sectors of fixed topological 
index with the answer obtained by summing over the topological charge. Also the microscopic 
Dirac operator spectrum in the full finite-volume Yang-Mills theory is obtained in this way, by 
summing over all topological sectors with the appropriate weight. 
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1 Introduction 



Although the chiral SU l(./Vj)xSU ji(Nf) symmetry of QCD is not spontaneously broken in a finite 
volume, it has long been known that many low-energy aspects of finite- volume QCD nevertheless can 
be derived from an effective chiral Lagrangian M . This is because the dominating degrees of freedom 
can still be the "would-be" Nj — 1 pseudo-Goldstone bosons (which are simply called pions in what 
follows) of the presumed chiral symmetry breaking pattern of SUL(iVj) xSUR(iVy) — > SU(iVj) in the 
infinite volume limit. Clearly at least two scales are of relevance here: if the linear extent L of the 
four- volume V is much larger than the Compton wavelength l/m n of the pions, the effect of the finite 
volume is mild, and can be treated as a small correction to the low-energy results of the theory in 
infinite volume. What is at first sight much more surprising is the fact that even in the opposite limit, 
when the pion Compton wavelenght much exceeds the size of the box, the effective partition function 
still becomes of a very simple chiral Lagrangian type. The procedure to derive this effective theory 
is to treat the zero modes of the pion field as new collective coordinates in the low-energy expansion. 
When all the dust has settled, the resulting leading-order Lagrangian coincides with what one would 
have obtained by the naive procedure of simply discarding all derivative terms in the infinite-volume 
chiral Lagrangian 0. Independently of whether chiral symmetry is spontanenously broken or not, the 
pseudoscalar states we here generically denote by pions are guaranteed, by exact mass inequalities, to 
be the lightest flavor non-singlet mesons in QCD ||]. 

Naively one might think that an effective Lagrangian technique as described above can provide infor- 
mation only about the chiral order parameter X = (ipifj) itself, or just soft pion dynamics. However, 
as shown by Leutwyler and Smilga [||, low-energy constraints on finite- volume gauge theories can, 
remarkably, be used to provide very non-trivial information about quark and gluon degrees of freedom 
as well. Specifically, the large-volume identification of the original QCD partition function and the 
effective Lagrangian restricted to the region 

' « L « — , (1) 



carries a lot of exact information about the spectral correlators of the smallest eigenvalues of the Dirac 
operator. While this holds for the full partition functions involved, the scoop of ref. ||] was to note 
that the analysis simplifies considerably if one focuses on gauge field sectors of fixed topological charge 
v. The effective partition function then becomes analytically solvable. It is thus possible to do a direct 
analytical comparison between, on one side, the full QCD partition function, and, on the other side, 
the effective partition function. One of the most striking outcomes of this analysis was the derivation 
of exact spectral sum rules for the Dirac operator in QCD and QCD-like gauge theories ||, |J. 

Surprisingly, the analysis of the spectral properties of the Dirac operator in finite-volume gauge field 
theories at fixed topological charge || has turned out to have an alternative description in terms of 
Random Matrix Theory [||, |j. The root of this equivalence lies in the universality of the Random 
Matrix Theory predictions: there is a huge class of Random Matrix Theories that, in the appropri- 



ate scaling limit, yield exactly the same spectral correlation functions |7| |8|, 10, 11]. Remarkably, 
so-called chiral versions of the three classical Random Matrix Theory ensembles have been shown 
by Verbaarschot to provide precisely the needed three universality classes of chiral symmetry break- 
ing in gauge theories with different gauge groups and different color representations of the fermions 
fl^| . While there are obvious resemblances here to the classification of critical behavior in terms of 
universality classes in statistical mechanics, the analogy should not be pushed too far. There is no 
renormalization group at work here, and the way QCD fits into one of these three universality classes 



1 



is by means of algebra (or, at a deeper level, symmetry) rather than dynamics. Indeed, within the last 
year it has become well understood how all the universal predictions derived from Random Matrix 
Theory can be rephrased in terms of the effective QCD partition function, suitably extended with 
additional quark species O]. Even universality of the results can easily be understood directly from 



the universality of the Random Matrix Theory partition functions fll4| and the equivalence between 
the effective chiral Lagrangian and the specific Random Matrix Theory with gaussian distribution, 
which was established in ||. By means of partially quenched chiral Lagrangians, it is now also possible 
to derive directly from the effective QCD partition function the microscopic spectral density of the 



Dirac operator, without any bypass through Random Matrix Theory at all [15|. All these more recent 
results permit us to say much more about the microscopic Dirac operator spectrum in finite-volume 
gauge theories. It is therefore appropriate to reconsider the analysis of ref. Q in this new light. 

It is thus partly by historical coincidence that almost all analytical predictions about the microscopic 
Dirac operator spectrum and the mass-dependent chiral condensate have been derived in gauge sectors 
of fixed winding number v. From the point of view of the "real" theory, where one sums over all 
topological sectors, this is analogous to knowing only the Fourier modes of a periodic function. While 
it may be interesting to study just these Fourier modes (and in the case of QCD much can indeed be 
learned from them), we would still rather like to know the full function. For QCD this means that we 
should know the predictions for the full theory at any given value of the vacuum angle 9, or, at least, 
since there are very stringent experimental bounds on the magnitude of this CP-violating #-angle, for 
the specific value of 6 = 0. 

The problem of summing over all winding numbers becomes particularly important when we consider 
tests of the analytical predictions for QCD sectors of fixed topological charge v by means of lattice 
gauge theory. Because lattice gauge configurations are not smooth on the scale of the lattice spacing, 
the very definition of gauge field topology is from the beginning an ambiguous concept. With so-called 
staggered fermions there is no help from index theorems (which otherwise can serve as an alternative 
definition of topological sectors) because there are no guaranteed exact zero modes away from the 
continuum limit In fact, lattice gauge theory data for QCD with staggered fermions and no 

restrictions on lattice field topology show remarkably good agreement with the analytical predictions 
for just the sector of zero winding number, v = 0. This phenomenon was first observed in comparisons 
of lattice data with analytical formulas for the mass-dependent chiral condensate (T^Jf], and later 
precisely the same effect was seen in lattice simulations of the microscopic Dirac operator spectrum 



with staggered fermions [18, [19| . One can argue that just because staggered fermions have no exact 
zero modes, this explains the puzzle. Are present-day lattice gauge theory simulations with staggered 
fermions thus only probing the sector of zero topological charge?Q Despite the problem mentioned 
above, lattice gauge field configurations generated by conventional Monte Carlo algorithms cannot 
all be classified as topologically trivial. They contain a distribution of configurations that smoothly 
fall into the precise topological sectors of the continuum theory as the lattice cut-off is decreased. 
Of course, in a given Monte Carlo sampling the sector corresponding to v = will be dominant. 
To understand the significance of just these configurations we need to know not only the analytical 
predictions for given sectors of fixed topological index, but also the weight with which these different 
topological sectors are contributing to full QCD. Very recently it has become feasible to address 
this in lattice gauge theory [2^j . But in the finite- volume region ([l]), these questions can also 
be answered analytically. For QCD-like theories the relevant partition functions at fixed topological 



1 Based on lattice data with dynamical quark masses so large on the relevant scale that the simulations effectively 
were quenched. 

2 In this connection, see the recent careful studies by the SCRI and Columbia groups [BOl. 
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charge v are known in closed analytical form [24]. It is simple to show that these fixed-f partition 



functions precisely give the non-trivial weights with which each fixed-f observable contributes to the 
full answer. Using this fact, we assess in section 2 the importance of the v = sector in bulding up the 
mass-dependent chiral condensate of the full theory. Section 3 is devoted to the microscopic spectral 
density of the Dirac operator, and we carry out the analogous procedure of summing over topological 
charges. Interestingly, one of the simple conclusions is that in the massless case the microscopic 
spectral density of the Dirac operator in the full theory simply equals that of the v = sector. This 
holds for all Nf > 2. The case Nf = 1 is - as expected - special, and there are then additional 
zero-mode contributions to the full microscopic spectral density from just the two sectors of v = +1 
and v = —1. In the massless limit just these two zero-mode contributions to the spectral density 
provide the required non- vanishing condensate for Nf = 1, even in this finite- volume situation. By 
numerically performing the sum over topological charges, we show how the full microscopic spectral 
density changes away from the v = prediction as the mass is increased for different values of Nf. 
We point out a difficulty with taking the quenched limit, and briefly discuss the more exotic chiral 
symmetry breaking cases, which correspond to different symmetries of the Dirac operator. All of the 
predictions described here can be tested in lattice gauge theory simulations. In fact, one of our main 
points is precisely that one can simultaneously test both the detailed analytical predictions in fixed-f 
sectors, and the predicted weight with which these sectors build up the full answer (which, in principle, 
should be much easier to measure). Section 5 contains our conclusions. 

2 Summing over Winding Numbers 



We begin with a few basic definitions. Consider QCD in Euclidean space on a finite 4-dimensional 
torus. Its four-volume can be taken to be just L 4 , or L 3 /T. The former case leads in the L —> oo 
limit to ordinary zero-temperature QCD, while in the latter case, if we take T fixed, the L —* oo limit 
gives Euclidean QCD at (small^) temperature T. For our purposes there is no need to distinguish 
between the two cases, as the infinite (three-)volume limit of the chiral condensate just must be chosen 
appropriately. The four-volume is in both cases denoted by V. As discussed at length in e.g. refs. 
II p5|] , gauge fields on such a four-torus can be described by a gauge potential A^x) on the full 
Euclidean space, plus a transition function Q(x), which describes the transition from one periodic cell 
to another (possibly non-trivial because A^(x) need only be periodic modulo gauge transformations). 
On the lattice one often prefers to work with strictly periodic gauge link variables (and hence trivial 
transition function), but by a gauge transformation these links can be made equivalent to a set of 
link variables with a non-trivial transition function. In the continuum such a gauge transformation 
will, however, necessarily be singular. Non-trivial winding numbers arise precisely from non-trivial 



transition functions [25]. 



For an SU(iV c ) gauge theory with Nf fermions in the fundamental representation, antiperiodic bound- 
ary conditions on the four-torus for these fermions the topological charge 



v - 

32v 



„ / d xTr [Fui/ Fpa\t uv per (2) 

In J 

runs over integer values. With a 0-term in the action, the partition function thus reads 

oo . 

Z(6) = Y, eiV(> [dA] u Y[det[0 - m f }e SYM[A] , (3) 



3 The temperature must be bounded by m T < T < Aqcd, in complete analogy with the condition AW) 
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after integrating out the fermions. Alternatively, 

oo 

Z(9) = ^ d Z v . ( 4 ) 



Z - = 7T- r dQe~ lve Z(d) . (5) 



where Z v is like a Fourier transform of the full partition function: 

1 

2vr Jo 

Consider first the case Nf > 2, and assume that chiral symmetry in the infinite- volume limit breaks 
spontaneously according to SlJ^iVy) xSUr(Aj) — ► SU(A r j). As shown in ref Q, the effective partition 
function then takes on a very simple form in the finite- volume region (|l]): 



Z{6) = [ dUexp \vZRe[e id/N fTrMU^}} , (6) 



>SU{Nf) 

where yVf is the quark mass matrix^, and E is the infinite-volume chiral condensate (for the massless 
theory with 9 = 0). One of the main lessons to learn from the analysis of ref. || is that these 
finite-volume effective partition functions for Nf > 2 in surprisingly many ways resemble the finite- 
volume partition function with Nf = 1. In fact, although the case Nf = 1 is radically different in that 
the axial U(l) symmetry is broken by the anomaly rather than spontaneously, the effective partition 
function turns out be given by what one naively could have extrapolated from (|6|), namely the simple 
exponential || 

Z{9) = exp rVERe[e w m]l . (7) 

In all cases the effective partition functions depend only on the combination m = EFwij. Keeping 
\ii fixed as the four-volume V is taken to infinity thus entails taking the chiral limit in a correlated 
manner. Moreover, the bound (jl]) must be satisfied throughout. In statistical mechanics such a 
procedure is known as finite-size scaling, and the \ii are finite-size scaling variables. 



The Fourier coefficients Z v have been computed analytically for all Nf ||, 24]. For Nf > 2 they can 
conveniently be written 

*,(W) = "^f 1 , W 
where the Nf x Nf matrix A({fj,i}) is {I n (x) is the modified Bessel function of order n) 

AM) = rf'^v+j-iM » ( 9 ) 

and 

AO. 2 ) = no* 2 -/*?) ( io ) 

i<j 

is the Vandermonde determinant. With an obvious interpretation of the formula (@), it even includes 
the case Nf = 1, for which || 

zm = iM . (ii) 

One notices that the general expression (||) is symmetric in v, 

MM) = z-v{M) . (12) 

as follows also directly from the original definition, eqs. (|) and §). 

4 From now on we always take this mass matrix to be real and diagonal: M = diag(mi, ...,mw,). 
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2.1 The chiral condensate for Nf — 1 



Having already observed that in the limit ([l]) the relevant scaling variables are /ij = T,Vm,i, we from 
now on try to eliminate as many factors of V and £ as possible. We hence measure all mass-dependent 
chiral condensates £(#; {//i}) in units of £, and effectively work with a unit four-volume. The mass- 
dependent chiral condensate for Nf = 1 is then simply defined by 

E(0;/i) = ^lnZ(0;v) . (13) 

In this Nf = 1 case the result is trivial and in fact //-independent: £(#; /i) = cos(#), as follows from 
eq. (]TT|) . 5 In any case, we shall mainly be concerned with the case 9 = 0, for which, in our units, 
£(0;/x) = l. 

Let us now compare with what we would find if we were to restrict ourselves to gauge field averages 
over fixed topological charge only. As already observed by Leutwyler and Smilga ||, this restriction 
is quite drastic. By definition, 

E„( M ) = ^-ln^Gu) , (14) 



which gives 



• ;/-!- ! w ' i ' ' i'\r i ■ (15) 

w) L /" J 

For v positive, the last term can easily be traced to the \v\ zero modes of the Dirac operator Z/*: In 
terms of the non-zero Dirac eigenvalues A n the original partition function (||) reads 

u=-oo f f n 

and the last term in ( |l5| ) precisely matches the term due to the prefactor m)"*. However, when v 
is negative this interpretation is not correct: in that case the sign of the second term in (O) is just 
the opposite of what is needed. There is no mystery here, but only a slight complication which 
we shall return to in section 3.1. Indeed, from the very definition (|14|) and the property (12), the 



mass-dependent chiral condensate (|lq ) is seen to be symmetric in v. S^(/u) = S_^(/i). Therefore the 
obvious ~ I/fx singularity in eq. (^) for v > is actually matched by a singularity of exactly the 
same strength for v < 0. 

We start by considering the case 6 = 0. In fig. 1 we show E^(/t) for a few small values of u, and 
compare each of them with the full answer £(/t). Contrary to the full answer which in this finite- 
volume limit is just the trivial constant E = 1, the chiral condensates in individual fixed-z/ sectors 
have non-trivial mass dependences. In this case we certainly cannot conclude that the v = sector 
adequately reproduces the full answer. But fixed- chiral condensates also show another disturbing 
feature, which we already mentioned above: for v ^ they diverge in the limit \i — ► 0. Actually, the 
second term in eq. ( |l5| ) by itself does not contribute to the full sum over topological sectors. To see 
this, we first derive a simple formula that tells us how to relate the full chiral condensate £(#; //) to 
the chiral condensates E„(/z) in sectors of fixed topological index: 

E(0; M ) = ^-hxZ(e ]f i) 



The ^-dependence of the condensate is not interesting; for 8 ^ it simply shifts from (tpip) to {ip^ tp) 
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We have here defined an "average over topology" by 

1 oo 
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Figure 1: The mass-dependent chiral condensate for Nf ~ 1. Shown are the results for a few fixed 

values of v. Even for 9 = they are quite different from the full answer (the constant 1 in that case). 



The relation (|17| ) tells us with which weight the chiral condensates in fixed topological sectors con- 
tribute to the full answer: Each S^(/u) comes with weight factor e lv9 Z v {fj). In fig. 2 we show how 
the full answer, the constant £ = 1 for 6 = 0, is built up as we sum over topological charge v from a 
given —\v\max to +\v\ max in eq. (|l7|). Perhaps surprisingly, we find that the v = contribution alone 
is not very significant, never exceeding more than around 20% of the full answer. In the shown range 
of fi we also see that by summing from, say, v = —8 to v = +8 we essentially obtain the complete S. 
The divergences that appear in £^(//) for v ^ are obviously completely removed when one performs 
the required average over topology. In the present case of 9 = this can easily be seen analytically. 
According to eqs. (15) and (17), and using Z(8 = 0;ii) = exp[/i] 

£(0 = 0;//) = e 



we here have 



oo 



' = — oo 
oo 

E 

'= — oo 



(19) 
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Figure 2: The full mass-dependent chiral condensate E(6* = 0,/i) for Nf — 1. The exact answer equals 1. 
Shown are the results of summing over the first few values of v. The v = sector is not particularly significant. 



The second term in the first line actually does not contribute at all: This is obvious for v = 0, and for 
v ^ the contribution from +v precisely cancels the contribution from —v. Next, substituting the 
Bessel function identity 

oo 

Y e iv6 I u {ii) = e" cos W , (20) 

v=— oo 

we of course recover the correct result £(0 = 0;/i) = 1 from eq. (|l9|). 

For a non-zero vacuum angle 9 the ^//U-term in £ v (/i/) does not disappear as simply when we sum 
over different topological sectors. Again this Nf = 1 case is very illustrative. Now, using some simple 
Bessel function identities, we have 

oo r 

V 



£(0;/x) = e^ cos(e) Y ^ 



e 



v=— oo 
oo 



A 1 



2 

cos(0) . (21) 



where in the last line we have agin used eq. (^0|). In this case the second term in (p^|) is certainly 



required to obtain the right answer from the sum (17). More generally, we also learn from eq. (17) 
that the fi — ► singularities in for v ^ are killed in the sum over topological sectors because 

these divergences are weighted by e lu6 'Z u (n). In general, Z u ([i) ~ n N f\ v \ for small /x. This follows both 
explicitly from the representation (|6|), or, in the original QCD language, from the Dirac determinant 
and the index theorem. Thus only for Nf = 1, and in that case only for the sectors with u = ±1, 
can these singularities in £ y (/x) give any non-zero contribution to the massless chiral condensate. 
In fact, one easily checks that in this case just these two contributions precisely yield the finite value 
H^O) = cos(0), as required H. 
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2.2 More flavors: N f > 2. 

It is convenient to first restrict ourselves to the case of equal masses, \ii = fi, Vi, so that the chiral 
condensates of all Nf species become equal. Then 



S(M = ±-^\nZ{e-ii) = ((£,(0;//))) , 



(22) 



with T, u (6; fj,) defined analogously. Consider first the case Nf = 2. The full partition function is known 
analytically ||: 

= f =, 1^2(1 + 008(6))^) . (23) 



1 . 4 
1.2 

1 

. 6 
0.4 
0.2 





1 '. 1 
I 


1 1 1 

v=0 


\ 


v=l 


\ 

\ 


v=2 - 


\ 


E 


\ 




\ 




\ 

s 




1 


i i i 



10 



Figure 3: The chiral condensate Y> v (p) for Nf = 2 and equal masses [i. There is no qualtitative difference 
with the picture from the Nt = 1 theory. For comparison we have also displayed the full chiral condensate 
for 9 = 0. 



Let us again first restrict ourselves to the case 6 = 0. The full mass-dependent chiral condensate then 
becomes 

J (2/x)+J 2 (2/x) 1 
S(0 = O '^ )= 2/ 1 (2 / x) 2V (24) 

which, despite appearances, is not singular at fi = 0. The partition function (||) reads 

Z„(/z) = I v (li? - Iu+lWv-M , (25) 

which gives 

v flA 1 f Iv(y)Iy+l(n) ~ J„-i(//)/„+ 2 (/i) , ^1 , 9R s 

" 2\ W-WMUM + M ' (26) 

We note again the presence of an explicit v / /i-term from the zero modes for v > (in which case the 
first term in (^) indeed remains finite in the limit ji — > 0). This term will always be present in E I/ (/x), 
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for any Nf. It reads in general u/(Nffj,), and thus contributes 



1 1 



£ uZ v (B = 0;n) 



to £(0 = 0; fx). Using the symmetry (|12|), this is seen to vanish: At 9 = this contribution to £(0; /i) 
always sums up to zero, for any Nf. When 9^0 this is no longer true: 

J r~ u= — 00 j r~ i/= — oo 



1 d 



N f ii 39 

In particular, for the case at hand, Nf = 2, we find that this equals 

ism(9) 



(27) 



2/z(l + eos(0)) 



/ 2 (/i(2(l + cos(0)) 



1/2N 



(which indeed vanishes at = 0). Moreover, we notice that this term is always strictly imaginary. 
Because the full sum is guaranteed to be real, this means that in some sense even for 0^0 the 
contribution of this term to T, u (9; /i) is only apparent; it only serves as to remove a corresponding 
imaginary part from the rest of the sum, leaving a remaining real answer. 
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Figure 4: The full chiral condensate S(/i) for Nf — 2, 6 — 0, and equal masses /x. We again compare 
different approximations, by summing only over the first few topoloigcal sectors. The v = sector is again not 
particularly significant. 



We now check to which degree of accuracy the chiral condensate £(#;//) for Nf = 2 and 9 = (i.e. 
eq. (|2"4|)) coincides with just the chiral condensate in fixed- sectors (|26|). This is shown in fig. 3 for 
a few small values of v. The divergent zero- mode contributions of course make the condensates for 
v 7^ very different from the complete answer. And actually in this Nf = 2 case the chiral condensate 
in just the v = sector at least has the qualitatively right behavior. But again, to see to what extent 
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the different topological sectors contribute to the full answer, we must add up the different sectors 



with the right weight. In fig. 4 we show how the full chiral condensate (24) is built up as we sum over 
topological charge v from given —\v\ max to +Hmax in eq. (gg). As in the Nf = 1 case, we find that 
the v = contribution alone is not particularly significant at all. It is obviously a generic feature, 
which is not specific to the very special circumstances of Nf = 1. The cases with Nf > 3 are indeed 
not qualitatively different [26|, and we do not display any of them here. 



3 The Microscopic Spectral Density 

Eigenvalues A n of the Dirac operator are defined by 0<p n = X n cj) n , and the spectral density is then 

= (E^ A " A «)) > ( 28 ) 

n 

where the average includes the sum over topological sectors. Because the weighting is done with the 
fermion determinant, the density depends on the masses m too. The double-microscopic limit is taken 
by keeping all pi = VErrii and £ = V£A fixed as V — > oo. The definition (|28| ) of course extends to 
averages in fixed-^ sectors as well: 

pM(A;{/Zi}) ^ { J25(X-X n )) u (29) 

n 

We now go to the double- microscopic limit. As for any other gauge-field expectation value in fixed- f 
sectors, also the microscopic spectral density p$ (Cj{Mi}) can be summed to give the microscopic 
spectral density of the full theory: 

oo 

ps(C-A{ih}) = £(0; {Mr 1 E ^MMWfoM) ■ (30) 

u=— oo 

It is useful to separate out the zero-modes explicitly: 

oo 

ps(C;0,M) = mMr 1 E ^MM)[p ( s\<:-,{fM}) + W\s{.0[ 

u=—oo 

= psiC-AM) + ((\u\))6(0 , (31) 
with an obvious definition of ps(C) ®i {pi})- 

It has been shown ]l3|, 15] that the microscopic spectral density (p^) without zero-mode contributions 



can be expressed directly in terms of the effective partition function, 

#(C;{M) = ~(-iriCI II(C 2 + m^ - J.) C ' • ( 32 ) 

f=i Zu ({in}) 

where we for clarity have explicitly indicated the number of fermion species that are involved on the 
right hand side. Actually, from the connection to Random Matrix Theory it is also known that all 
microscopic spectral correlators Ps ■ ■ ■ , ( n ) can be derived from just one "master formula" 

m, 

1 N f , Z {Nf+2) (\u \ iC iC) 

/=i z » \\m) 
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All spectral correlation functions follow from this single function: 

P { s\&, ■■■,&,{*}) = detKP(Ca,Cb\{lM}) , (34) 

a,b 

where the determinant on the right hand side is taken in terms of the n x n matrix in the arguments 

(a, Cb- 

The full microscopic spectral density without zero-mode contributions, ps(C'i ®i {pi})> takes a particu- 
larly simple form in the massless limit. Using Zq(Q) = 1 and Z u (0) = for all v ^ 0, we immediately 
get 

ps(C-Mih = o}) = P ( s = °\C;{Pi = o}) • (35) 

That is, in the massless limit the full microscopic spectral density without zero modes simply coincides 
with the corresponding density of the v = sector. Note in particular that in this massless limit the 
full microscopic spectral density is independent of the vacuum angle 9. This is as expected: in QCD 
with at least one massless flavor any #-angle can be removed by a chiral rotation. 

The result ( |35| ) may seem to contradict the fact that in the Nf = 1 theory the chiral condensate is 
= 0) = cos(#) even though the massless limit is taken at finite volume. Indeed, even at finite 
volume there is a kind of Banks-Casher relation for the microscopic spectral density too, 

lim(2 M ) rdt PS ff'f = 7r Ps (0;6,0) = Trp^OjO) = 0. (36) 

This microscopic spectral density alone can therefore never give rise to a condensate in the massless 
finite-volume theory. But the full condensate potentially also receives contributions from the zero 
modes, even at zero mass. This part equals (taking Nf equal masses for simplicity) 

oo 

(Nfpr'dH)) = [NfpZie;^})}- 1 £ e**MM)W\ 

u=—oo 

oo 

= 2[iV // iZ(^{/i i })]- 1 ^cos(^)^({/x i }) l / . (37) 

v=\ 

In the limit p — > this expression vanishes for all Nf > 2. In view of ( |36"D this just means that for 
Nf > 2 there are no chiral condensates in the massless limit, in agreement with the fact that there 
cannot be spontaneous symmetry breaking in these finite-volume theories. For Nf = 1 and 9 = the 
sum in (|37| ) can be performed analytically, giving || 

p- 1 m))=e-»[I Q {p) + h{p)] , (38) 

whose limit at p = indeed equals unity, as required in order to obtain S = 1 in that case. In general, 
for 9 we note again that because of the behavior I n (x) ~ x' n ' as x — > 0, only the v = ±1 terms 
contribute to the sum when p — > 0, yielding precisely the required — > cos(#) in the limit 
H-> 0. 

Inserting the formula (^) into the definition of ps(C] 9, {pi}) in eq. ([H]) gives an interesting compact 
expression for the full microscopic spectral density: 

oo 

PsiC-AM) = zW(9-{p t })- 1 £ e*"MM)p%\C;M) 

v=—oo 
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-1 N f oo 



/=i 

AT/ 

icill(c 2 + ^) 

/=1 



Z( N f+ 2 \e + Tr;{^},i(,i() 
Zi N f){9-{^}) 



(39) 



While this formula strongly resembles the one for fixed we note the effect of the (— l)^-factor in the 
individual p^(A; {/UjD-contributions: The partition function with two additional fermion species (in 
the numerator) is to be evaluated at a shifted vacuum angle: 6 + it. Thus even if we insist on getting 
the physical microscopic spectral density with 9 = (or at least exceedingly small), we need to know 
also the effective QCD partition function with two more flavors, and with a non-vanishing vacuum 
angle. 

One can perform a similar summation over topological charge of all microscopic spectral correlators 



Ps (Cl) ■ • ■ i Cn! {a*i})- However, the representation ( J33| ) does not give a a convenient compact expres- 
sion, due to the fact that the determinant must be taken inside the summation over v. Fortunately, 
there exists a different expression for Pg\(,\, . . . ,£ n ; {pi}) (see the third paper of ref. [13|), which, 
although more unwieldy in general when v is fixed, is more convenient in this context: 



N f 



PsH(i, {M) 



c^u kn«?+* II(cJ- 

i \ f=l ) j<l 

x ^ +2fe) ({M,Ki},---,Kfc}) 



2\2 



(40) 



Summing this expression over v will yield the microscopic spectral correlators of the full theory, now 
expressed in terms of the ratio of a partition function of Nf + 2k species and the usual partition 
function corresponding to the actual number of fermions Nf . 
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Figure 5: The full microscopic spectral density of the Nf 
with the v = result. 



1 theory with 6 = 0. The massless case coincides 
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The full effective partition function (^|) is unfortunately not known in closed analytical form for Nf > 3. 
So to evaluate the right hand side of eq. (^) it is most convenient to simply perform the sum over 
topological charges numerically. We know already the full answer in the massless case (|35|), and it is 
thus interesting to trace the change in ps(\;8,{fii}) as the mass is increased. In fig. 5 we show the 
full microscopic spectral density for the Nf = 1 theory with 9 = 0. The /x = curve indeed coincides 
exactly with the massless prediction of just the v = sector H 



("=0)/* 



=o) 



id ji(c) 2 -j 2 (c)«/o(c) 



(41) 



As the mass increases, sectors with higher values of \v\ begin to contribute. Because the oscillations in 
Ps (C; p) are shifted to the right as \u\ increases, there is a certain amount of destructive interference, 
and the full microscopic spectral density becomes smoother as /U is increased. The limit /U — > oo, which 
naively should correspond to the quenched theory, is ill-defined, as will be discussed below. 

. 5 
.45 
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^1=0,^2=0 
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Figure 6: The full microscopic spectral density of the Nf = 2 theory with 9 = 0. When both fermion masses 
are set to zero we recover the v = result. When one of the masses is taken to infinity, we recover the full 
microscopic spectral density of the Nf = 1, by decoupling. 



The case Nf = 2 is more interesting. Because of the simple identification (p5|), the \x\ = ^2 = limit 
simply reproduces the massless Nf = 2 density of the just the v = sector, 



0=0) /> 



= /X2 = 0) 



2ICI [MO 2 



(42) 



0, 



In another extreme limit, when, say, [12 — * 00, this second fermion decouples. If we still take /ji 
we thus recover the simple massless prediction of just the v = sector of the Nf = 1 theory, eq. (41). 



For intermediate masses, there is again some amount of destructive interference that smoothes out the 
microscopic spectral density. Shown in fig. 6 are the two limits mentioned above, and one example of 
the prediction for intermediate masses. 
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3.1 Careful with the sign of v 



We have already noticed that the symmetry (12) under v — ► —v also implies S v ({//j}) = S_ I/ ({//j}). 
But this holds only if we consider all contributions, including the v/p-piece. For example, for Nf = 1 
(and 8 = 0) the two terms on the right hand side of eq. (|l^) are not separately invariant under 
v — > —v, but the sum is. This has consequences for our definition of the microscopic spectral density 
Pg(\, {/Xj}), which does not include the contribution from zero modes. 

It follows from eq. (12) and the representation (32) that also 

pPMlH}) = P^MlH}). (43) 



Clearly also the full microscopic spectral density (|29j) is thus symmetric under v — ► —v. Let us for 
illustration consider one of the "massive spectral sum rules" for Nf = 1, which was derived in ref. |2(| 
for the case v > 0. First, by definition, 



/■OO 

Sj/(/x) = 2/i / dX 
Jo 



#(A;/*) 



Comparing this with eq. (|l5|), we see that there are actually two different massive sum rules, depending 
on whether v is positive or negative: 

u>0: HdX^^ - ^ 



A 2 + /j 2 2nl v (n) ' 

,<0: r d X P t^4 = ^44 + 4- (45) 

Jo A 2 + /i 2 2fjJ v (jj) fi 2 K ' 

It is straightforward to check that in the limit \i — ► these two relations precisely lead to two massless 
spectral sum rules that can be combined into (the sum running over strictly positive eigenvalues only): 

(?'£). -Jwnj- w 

in complete agreement with ref. These considerations of course trivially generalize to any Nf, 
where only the right hand side of eq. (051) becomes increasingly more complicated 



The Nf = 1 case that we have used as an example here is particularly simple also when it comes to 
an analytical check on our formula for the full microscopic spectral density (p9|). In that case one can 
explicitly confirm the relation 

= 2 M /V P ffi C, £° +-<<M)) (47) 

since all the involved integrations can be performed analytically. (The left hand side is of course just 
the trivial constant cos(#) in this case, and the full spectral density ps(0] C> A 4 ) I s indeed just a function 
that achieves this ^-independence of the right hand side as well). 
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3.2 Trouble with the quenched limit 



The quenched limit of Nf — > is not easily taken. In gauge-field sectors of fixed topological charge v, 
the quenched mass-dependent chiral condensate reads |l^, [l^] 

£„(//) = n (I v {n)K v (n) + I u+1 (n)K u ^)) + - , (48) 

A 1 

where K n (x) is the modified Bessel function. This quenched chiral condensate also satisfies T, v (jjl) = 
S— It shares the ~ l//i divergence with the unquenched condensates when ^ ^ 0. In the latter 
case we have already seen that this divergence is harmless in the sense that it is killed once we sum 
over topological charge. Here this is no longer true. If we simply assume that Z v (n) = 1 for all v in 
the quenched limit, it is not even possible to sum over topological charge. There is now no suppression 
of higher-i^ contributions, and the sum in 

1 oo 

EOM = zTKTTK E ^E,( M ) (49) 

does not converge. If one tries to "regularize" the sum by introducing some maximum topological 
charge I 

v\max &nd defining 

Z(0;v) reg = ]T e ive E v (fx) , (50) 

eq. (^) still has no well-defined limit as |z/| m<M; — > oo. This failure of convergence is not related to 
the fj, — * singularities in fixed-i^ sectors; it arises here from the infinite summation over topological 
charge, if, in the quenched approximation, we assume that each topological sector enters with equal 
weight Z y {^i) = 1. It is clearly invalid to make this assumption. In fact, in the pure gauge theory 
the topological susceptibility \t = ((z^ 2 )) ( m our units) is a well-defined and finite number, and the 



distribution of topological charge is believed to be Gaussian in the infinite- volume limit [27]: 



Z, 



V 



e~ v V 2 ** . (51) 



Z 2ir X t 

It is tempting to simply insert this ansatz for Z v into eq. (^) instead of unity, but it is not quite 



consistent. The whole framework here is based on chiral Ward identities, of which [28| 



N f ((u 2 }) = n (52) 

indeed is a direct consequence of these finite- volume partition functions Q . Taking the quenched limit 
detaches us from this tight constraint, and a different analysis is required to establish the behavior 
of Z v as a function of i>. In any case, the /i — > singularities are no longer removed by the sum 
over topological charge, because the distribution of Z v will be //-independent. This ~ 1/fi singularity 
in the (full) quenched chiral condensate has recently been confirmed by direct lattice measurements 
using both domain-wall and overlap fermions (2C 



With these observations, we should certainly expect difficulties with the quenched microscopic spectral 
density. The quenched microscopic density in a gauge field sector of fixed v reads || [lTj 

?f{Q = [UC? - <WCR-i(C)] • (53) 
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According to eq. (31), and if we again take Z u (fi) to be a ^-independent constant (say, unity), the 
sum in 

oo 



(54) 



again requires regularization. Introducing a maximum topological charge | 

v\max — * oo as above, we 

actually find a remarkable regularity. Let us for convenience restrict ourselves to 9 = 0. Shown in fig. 
7 is a series of successive summations with increasing 1^1 max* The resulting function is approximately 
iineaTrL until, at a point that shifts to infinity with increasing 1^1 max? it turns over into the constant 
value 1 /tt. The slope of the linear part thus depends on the choice of | 

v\maxi and there is no convergence 

of the sum. As discussed above, the approximation of taking all Z v to unity is just drastically wrong. 




25 30 



Figure 7: Summing over topological sectors with incorrect weight: if in the quenched theory one simply 
assumes that Z v = 1 the resulting microscopic spectral density becomes linear up to the value 1/tt, where it 
turns flat. This point where the linear part is replaced by a flat part moves to infinity with increasing |^| maa; , 
and thus has no limit. 



Another way of reaching the quenched limit is by decoupling all Nf fermions in any given theory by 
sending all masses to infinity. For simplicity taking Nf = 1, and sending the mass \x to infinity, we 
obtain precisely the same kind of curve as obtained above in the naive way of just summing over each 
fixed-z^ sector with equal weight. The cut-off which controls the point where ps(C) turns from linear 
to flat behavior is now played by this mass [i. This suggests that decoupling all fermions by taking 
all their rescaled masses /Xj to infinity is not providing a physical limit of the full theory. Indeed, 
one easily checks that in this extreme limit all partition functions Z v become independent of v, thus 
explaining why the two different ways of trying to obtain the quenched microscopic spectral density 
described here fail in precisely the same manner. 

How does this result fit in with the general formula (^)? Taking the Nf — > limit in that equation 
6 I thank S. Nishigaki for first pointing this out to me. 
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gives 



1 Z( 2 Xe + n;i(,iQ 

ps(x-,e) = -\c\ 



(55) 



For the quenched microscopic density we thus need the effective partition function for Nf 
fortunately is known in analytical form (see eq. (23 
indeed find from eq. (|23|) that 



2, which 

For simplicity taking the case of 6 = 0, we 



P s(x-,e=o) 



-ICI 



(56) 



2Z(°) 1 

The density consistently comes out linear, and again with an undetermined constant in front. To 
determine this constant we would normally impose the matching condition of ps(( —> oo) = 1/tt, an 
impossible requirement in this case. Again, the failure of these attempts at taking the quenched limit 
in this framework is perhaps not surprising. The presumed distribution ( |51~D concerns by construction 
only the pure gauge sector, for which no information can be gained via the (anomalous) chiral Ward 
identities on which the present framework is based. Moreover, in the pure gauge theory we can sum 
over gauge field configurations with also fractional winding numbers pa]. In the large- N c limit the 



topological susceptibility in the quenched theory is, via the Witten-Veneziano relation [29], related to 
the flavor singlet rj mass. In the quenched limit, where the anomaly plays no role, the rf should be 
treated as a (pesudo)-Goldstone boson too. Perhaps an analogue of the large- N c effective Lagrangian 
with built-in information about the topological susceptibility of the pure gauge theory (as required 
by the Witten-Veneziano relation, at least at large N c ) can be used to treat this case in a consistent 



manner. 




Figure 8: The partially quenched chiral condensate £(/x„,/i) in the Nf — 1 theory with 6 — 0. The 
"quenched" mass is indicated by /i v , while /i is the real physical mass, here taken to be unity for convenience. 
The sum over topological charges converges everywhere, except at the point \x v = 0. In this case the divergence 
at fi v — > can easily be worked out analytically (see eq. (|5^)), and for comparison we display also just the 
resulting term with the coefficient as given in eq. (KJ). 



There are no such difficulties with partial quenching. For example, the dependence of the partially 
quenched chiral condensate on the (quenched) "valence quark" mass \x v is computed in the full theory 
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by averaging with respect to the usual fixed-f partition function Z v . This sum converges for fj, v ^ 0. 
To be concrete, consider the analytical expression for the Nf = 1 case which was derived in the second 



ofref. gj: 



| 2 ^ K ^v) fiyIu(fiv)Iu+l(lJ-) ~ llI v (n)I v +l{Hv) | V ^ 

I v (fx) n% - /i 2 fa 

In fig. 8 we show this as a function of fi v for v = (which of course does not exhibit the l/fi v - 
singularity) , and fixed value of the physical "sea quark" mass fi. We also show the full answer, for the 
same mass fj,, obtained by summing over winding numbers according to 

oo 

= 2(9; fi)- 1 e^OOE^M) , (58) 

u=— oo 

where in this case Z(8\\i) = exp[/tcos0], and Z u (fi) = 

The partially quenched cases are particularly interesting because they provide a limit in which the 
~ 1//Jiv divergences of the quenched approximations are under full control from the analytical point 
of view. Indeed, the coefficient of the l/^-term in the partially quenched Nf = 1 theory can be 
computed analytically in this finite- volume limit; it equals 

((H)) = ^-"(loM+IiOi)). (59) 

This curve has also been included in fig. 8, and one clearly sees how accurately it describes the l/fi v 
singularity. 

3.3 Extensions 



The present results easily extend to the other two major classes of chiral symmetry breaking [12|. The 
analytically most simple case is that corresponding to gauge group SU(iV c > 3) and Nf fermions in 
the adjoint representation (which in Random Matrix Theory language corresponds to the chSE). For 



example, from the relations derived in the second of reference [13] we know that we in this case can 



express the microscopic spectral density of the Dirac operator in a sector of topological charge v by 

#'(C;W) - CC ^"''™ . (60) 

ffl ''({«}) 

Also this can be summed over topological charge v to get the full microscopic spectral density for that 



case. Using a similar relation for all higher /c-point functions [13 



k ( N f \ k 

P { s\ci,...,t k ;{»i}) = cf n cfn«?+M/) nici-cfi 4 

^ /+4fe) ({M;Ki},---,Kfc}) 



(61) 



also these fc-point functions can be averaged over topological charges to yield a ratio of two full partition 
functions, without restrictions to topology. In these relation each (imaginary) additional fermion mass 
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on the right hand side is four-fold degenerate. Unfortunately, the relevant fixed-f partition functions 
are not in general known in closed analytic form, and the proportionality constants in both of the 
above relations have therefore not yet been determined. Until this has been done, and in particular 
their ^-dependences have been found, we cannot explicitly identify the combination of finite-volume 
partition functions that will emerge from the sums. 



4 Conclusions 



We have here considered the effect of summing over all topological charges u, while remaining in the 
finite- volume scaling regime defined by eq. ([l]) . For the mass-dependent chiral condensate the effect of 
summing over topological sectors is quite drastic. For example, the condensate actually diverges like 
~ 1/fi in the /i — > limit if one restricts oneself to any sector of non- vanishing topological charge v, 
even though the full condensate vanishes as \x — > (except for the case Nf = 1). The microscopic Dirac 
operator spectrum itself is also substantially modified by the summation over topological charges. A 
simple compact relation gives the full microscopic spectral density in terms of a ratio of partition 
functions, one of which is evaluated at a shifted vacuum angle 6 — > 8 + tt, and with two additional 
fermion species. We have pointed out a difficulty with obtaining the quenched (Nf — > 0) limit in this 
simple way. The difficulty disappears if one focuses instead on a partially quenched limit, in which 
the appropriate weight factors in the sum over topological numbers again are known exactly in this 
finite-volume regime. We have also noted a drastic simplification in the massless case, where the full 
microscopic spectral density actually coincides with the one of the sector with v = 0. 

An interesting open question is whether one can establish a compact Random Matrix Theory formu- 
lation also for the full effective partition function Z(8; {/Zj}). Although we believe the answer to be 
negative, we are not aware of any strong argument to that effect. 



In future lattice studies it seems obvious to try to study not just the microscopic Dirac operator 
spectrum in sectors of fixed u, but also the associated distribution of configurations labelled by the 
topological index. Thanks to very recent developments this is now possible [21, 22]. This way one can 
test simultaneously the probability distribution of configurations with fixed topological index (which 
is known analytically in the regime (||) and given, for 6 = 0, simply by Z 1/ ({fii})/Z(9 = 0;{/ij})), 
and, simultaneously, the full microscopic Dirac operator spectrum. Such numerical simulations with 
dynamical fermions now seem feasible. 
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